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This paper reviews the recent developments on building nanoelectronics for our future informa¬ 
tion processing paradigm using multiferroic composites. With appropriate choice of materials, 
when a tiny voltage of few tens of millivolts is applied across a multiferroic composite, i.e. a piezo¬ 
electric layer stain-coupled with a magnetostrictive layer, the piezoelectric layer gets strained 
and the generated stress in the magnetostrictive layer switches the magnetization direction be¬ 
tween its two stable states. We particularly review the switching dynamics of magnetization 
and calculation of associated metrics like switching delay and energy dissipation. Such voltage- 
induced magnetization switching mechanism dissipates a minuscule amount of energy of only '^1 
attojoule in sub-nanosecond switching delay at room-temperature. The performance metrics for 
such non-volatile straintronic devices make them very attractive for building not only memory 
devices but also building logic, so that they can be deemed suitable for computational purposes. 
Hence, multiferroic straintronics has profound promise of contributing to Beyond Moore’s law 
technology, i.e. of being possible replacement of conventional charge-based electronics, which is 
reaching its performance limit specifically due to excessive energy dissipation. 
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1. Introduction 

The conventional charge-based electronics for more 
than past fifty years has a history of great succesJI^. 
However, the proven concept of enhancing the per¬ 
formance metrics by miniaturizatioil^ of devices is 
approaching its fundamental limiti^, while there are 
issues due to process variation, basically the exces¬ 
sive energy dissipation in the devices limits the fur¬ 
ther improvement of transistor-based electronicJ^. 
The Nanoelectronics Research Initiative (NRlj^ at 
United States says “Future generations of electron¬ 
ics will be based on new devices and circuit architec¬ 
tures, operating on physical principles that cannot 
be exploited by conventional transistors. NRI seeks 
the next device that will propel computing beyond 
the limitations of current technology.” The chal¬ 


lenge is to invent switching devices that dissipate 
miniscule amount of energy, e.g. ~1 attojoule while 
maintaining sub-nanosecond switching delay. The 
performance metric switching delay is particularly 
important because if transistors were to switch at 
low switching speed there would not be excessive en¬ 
ergy dissipation due to delay-energy trade-off. This 
is quite explicitly pointed out by NRI while seeking 
new switching devices over traditional transistors 
since several proposals could not meet the require¬ 
ment of switching speeds. 

Devices made of electron’s spin as a state vari¬ 
able are switched by flipping spins without moving 
any charge in space and causing current flow. This 
eliminates the ohmic loss associated with switch¬ 
ing, although some energy is still dissipated in flip¬ 
ping spins from one state to another overcoming the 
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switching barrier between the states. It is widely 
believed that using “spin” as state variable is ad¬ 
vantageous over the charge-based counterpart. Un¬ 
fortunately, however, this advantage will be squan¬ 
dered if the method adopted to switch the spin is so 
energy-inefficient that the energy dissipated in the 
switching circuit far exceeds the energy dissipated 
inside the switch. Regrettably, this is often the case, 
e.g. switching spins with a magnetic fielci^ ® or 
with spin-transfer-torque mechanisrrPl using an ex¬ 
ternal charge currenl!^. Hence, there is a need for 
inventing a switching mechanism to flip spins for 
the paradigm of spintronics to be technologically 
viable. 

According to Brown’s fundamental fine- 
ferromagnetic-particle theorjffSI, magnetic domain 
formation should be limited to very small dimen¬ 
sions (~100 nm) because of the competition be¬ 
tween the magnetostatic energy and the quantum- 
mechanical exchange energy, causing nanomagnets 
to behave like single giant spins. These “giant” 
spins can beat superparamagnetic limit at room- 
temperatur^m, which is crucial for general-purpose 
information processing. The minimum energy dissi¬ 
pated to switch such a single-domain nanomagnet 
(a collection of M spins) can be only ^kTln{l/p), 
where T is temperature and p is error probability, 
since the exchange interaction between spins makes 
M spins rotate together in unison like a giant clas¬ 
sical spirff^l ttSl contrary, the minimum 

energy dissipated to switch a charge-based device 
like a transistor would be ^NkTln{l/p), where N is 
the number of information carriers. This gives nano¬ 
magnets an inherent advantage over the traditional 
transistors with regards to energy dissipation. 

In multiferroics, different ferroic orders such 
as ferroelectric, ferromagnetic/ferrimagnetic, fer- 
roelastic etc. coexist. For our discussion, we will 
assume the coexistence of ferroelectric and ferro¬ 
magnetic orders to mean mutiferroism. Single-phase 
multiferroic materials are rare and moreover the 
magnetoelectric responses of those is very weak or 
occur only at low temperatures so their technologi¬ 
cal applications is not yet very promisingly. On the 
contrary, 2-phase multiferroic composites consist¬ 
ing of magnetostrictive layers strain-coupled with 
piezoelectric layer ^ ^ ^3 [IS j^ave the 

bottlenecks as of their single-phase counterparts. 
Thus multiferroic composites are more promising 


for technological applications. 




Magnet's plane: 4) = ±90° 


Fig. 1. A 2-phase multiferroic nanomagnet in the shape of 
an elliptical cylinder is stressed with an applied voltage via 
the dai coupling in the piezoelectric. The multiferroic is pre¬ 
vented from expanding or contracting along the in-plane hard 
axis (y-axis), so that a uniaxial stress is generated along the 
easy axis (a-axis). (Reprinted with permission from Ref. 19. 
Copyright 2012, AIR Publishing LLC.) 


The magnetization of the shape-anisotropic 
single-domain magnetostrictive nanomagnet in mul¬ 
tiferroic composites can be switched in less 
than 1 nanosecond while dissipating only ~1 
attojouljy [ISl. Hence such devices have emerged 
as potential candidates as storage and switching 
elements for our future non-volatile memory and 
logic systems. Particularly due to high switching 
speed while simultaneously being highly energy ef¬ 
ficient, this has lead to l ogic pr oposals incorpo¬ 
rating such systemJSni Hil |22] |23] magnetiza¬ 
tion of the nanomagnet has two stable states (mu¬ 
tually anti-parallel) along the easy axis encoding 
the binary bits 0 and 1. The magnetization can be 
switched from one stable state to the other when a 
tiny voltage of few tens of millivolts is applied across 
the piezoelectric layer while constraining it from ex¬ 
panding or contracting along its in-plane hard-axis 
(see Fig. [T]). The applied voltage produces a strain 
in the piezoelectric layer, which is then transferred 
to the magnetostrictive layer. This in turn generates 
a uniaxial stress in the magnetostrictive nanomag¬ 
net along its easy-axis and rotates the magnetiza¬ 
tion towards the in-plane hard axis as long as the 
product of the stress and the magnetostrictive coef¬ 
ficient is negative. It is assumed by convention that 
a tensile stress is positive and a compressive stress 
is negative. There have been experimental efforts to 
demonstrate such electric-field induced magnetiza- 
tion rnt.tinJ3Eliai2Zl|2a 


The rest of the paper is organized as follows. 
In Section |2l we describe the model for such mul¬ 
tiferroic straintronic devices. The Landau-Lifshitz- 
Gilbert (LLG) equation is solved to track the mag- 
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netization dynamics and to calculate the associated 
performance metrics like switching delay and energy 
dissipation. Also, we present the model for deter¬ 
mining magnetization dynamics in a circuit made 
of multiple multiferroic devices. Section [3] presents 
the simulation results for both a single memory de¬ 
vice and a circuit of multiple multiferroic devices. 
Finally, Section U] summarizes this review and pro¬ 
vides the outlook of the multiferroic straintronic de¬ 
vices on building nanoelectronics for our future in¬ 
formation processing paradigm. 

2. Model 

In this Section, we will hrst review the model of a 
single multiferroic device. The emphasis would be 
on magnetization dynamics of the magnetostrictive 
nanomagnet in a multiferroic composite bysolv- 
ing Landau-Lifshitz-Gilbert (LLG) equatioiPSl EQ] 
We will review the model incorporating room- 
temperature thermal fluctuations using stochastic 
LLG equatioEpD US. Then we describe how the 
same model can be used for unidirectional signal 
propagation in a chain of multiferroic devices^H 1221 
How we can build logic gates for computational pur¬ 
poses is described too in the end. 

2.1. Single Multiferroic Device 

We consider a single isolated nanomagnet in the 
shape of an elliptical cylinder with its elliptical cross 
section lying in the y-z plane; the major axis is 
aligned along the z-direction and minor axis along 
the y-direction (see Fig. [T]). The dimensions of the 
major axis, the minor axis, and the thickness are 
a, b, and I, respectively. So the magnet’s volume is 
H = (7r/4)a5L The z-axis is the easy axis, the y-axis 
is the in-plane hard axis and the x-axis is the out-of- 
plane hard axis. Since I <C b, the out-of-plane hard 
axis is much harder than the in-plane hard axis. Let 
9{t) be the polar angle and (j){t) the azimuthal angle 
of the magnetization vector in standard spherical 
coordinate system. Note that when (p = ±90°, the 
magnetization vector lies in the plane of the nano¬ 
magnet. Any deviation from (p = ±90° corresponds 
to out-of-plane excursion. 

We can write the total energy of the magne¬ 
tostrictive single-domain nanomagnet when it is 
subjected to uniaxial stress along the easy axis (ma¬ 
jor axis of the ellipse) as the sum of the uniax¬ 
ial shape anisotropy energy and the uniaxial stress 

lOOl 

anisotropy energy^. We assume that the magne¬ 


tostrictive layer is polycrystalline, so that we ignore 
the magnetocrystalline energy. The uniaxial shape 
anisotropy energy at an instant of time t is given 

b#2l 


EsHA{t) = EsHA{9{t),p{t)) 

= {go/2)M^QNd{e{t),cP{t)) ( 1 ) 

where Mg is the saturation magnetization and the 
demagnetization factor Nd{t) is expressed ad^ 


Nd{t) = Nd{0{t),(p{t)) = Nd-xxSin^9{t)cos^(p{t) 

± Nd-yysin^9{t) sin^(p{t) ± Nd-zzCos^9{t) (2) 

with Nd-ram being the m*^ (m=x,y,z) component 
of the demagnetization factoi®]. Note that these 
factors depend on the dimensions of the nanomag¬ 
net and not on the material properties. The dimen¬ 
sions of the nanomagnet is chosen as a = 100 nm, 
6 = 90 nm and I = 6 nm, which ensures that the 
nanomagnet has a single ferromagnetic domaiiff^l. 
These dimensions alongwith the material parameter 
saturation magnetization Mg determine the shape 
anisotropy energy barrier, which separates the two 
stable states of the nanomagnet. The in-plane en¬ 
ergy barrier Eg (p = ±90°, see Fig. [I]), which is the 
lowest difference between the shape anisotropy en¬ 
ergies when 9 = 90° and 0 = 0°, 180° determines 
the static error probability of spontaneous magne¬ 
tization reversal due to thermal fluctuations. Ac¬ 
cording to Boltzmann distribution, this probabil¬ 
ity is exp [—Eg/kT], This probability should be low 
enough for technological application purposes. With 
the dimensions and material chosen, Eg = 44 kT at 
room temperature, so that the static error proba¬ 
bility at room temperature is Note that the 

dynamic error probability on the other hand signi¬ 
fies the switching error probability when magneti¬ 
zation fails to flip from one state to another during 
switching events. Usually the dynamic error prob¬ 
ability can be much higher than that of static er¬ 
ror probability. So we must put emphasis on this 
dynamic error probability to meet the demand of 
technological viability. 

The stress anisotropy energy is given b; 


EsTA{t) = EsTA{9{t),a{t)) 

= —{2>/2)\sa{t)Q.cos^9{t) (3) 

where (3/2)A* is the magnetostriction coefficient of 
the magnetostrictive nanomagnet and apt) is the 
stress generated in it by an external voltage. A pos¬ 
itive Xgapt) product will favor alignment of the mag¬ 
netization along the major axis (z-axis), while a 
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negative \sO'{t) product will favor alignment along 
the minor axis (y-axis), because that will minimize 
EST Ait)- We will use the following convention that 
a compressive stress is negative and tensile stress 
is positive. Therefore, in a material like Terfenol-D 
that has positive A*, a compressive stress will favor 
alignment along the minor axis, and tensile along 
the major axis. The situation will be exactly oppo¬ 
site with nickel and cobalt that have negative A^. 

At any instant of time t, the total energy of the 
nanomagnet can be expressed aJl^ 

E{t) = Eie{t),(j){t),a{t)) 

= i?( 0 (t), + C'(t) ( 4 ) 

where 

Bit) = Bi(t)it),ait)) = Boicpit)) + 

^stress i(^it)) 
(5a) 

Boit) = Boicpit)) = iRo/2) M^n[Nd-xxCos‘^4>it) 

+ Nd-yySin^4>it) - Nd-zz] (5b) 

Egtressit) — B stress (fT(t)) = (3/2)AsO-(t)f] (5c) 

Cit) = Ciait)) = iRo/2)M^nNd-zz 

- (3/2)A,a(t)fi. (5d) 

The magnetization M(t) of the nanomagnet 
has a constant magnitude at any given temperature 
but a variable direction, so that we can represent it 
by the vector of unit norm nm(t) = lV[(t)/|lV[| = Or 
where Or is the unit vector in the radial direction in 
spherical coordinate system represented by (r,0,0). 
The other two unit vectors in the spherical coor¬ 
dinate system are denoted by eg and e<^ for 6 and 
(p rotations, respectively. The torque acting on the 
magnetization per unit volume due to shape and 
stress anisotropy iJ^^ 

TE(t) = -n^it) X VEieit),<Pit),<jit)) 

= —2Bi(j)it), ait))sin9it)cos6it) e^ 

- BQei4>it))sin0it) eg, (6) 

where 

Boeit) = Boei4>it)) 

= ( mo / 2 ) M‘^niNd-xx - Nd-yy)sini2(j)it)). 

(7) 

The effect of room-temperature random ther¬ 
mal fluctuations is incorporated via a random mag¬ 
netic field h(f), which is expressed as 

hit) = hxit)e„ + hyit)ey + hzit)e^, ( 8 ) 

where /ij(t) (i=x,y,z) are the three components of 
the random thermal field in Cartesian coordinates. 


We assume the properties of the random field h(f) 
as described in Ref. EH The random thermal field 
can be written aJ^ 


hiit) = 


2akT 

lyjMyAt 


<^{ 0 , 1 ) (i) ii = x,y,z) (9) 


where a is the dimensionless phenomenological 
Gilbert damping parameter, 7 = 2iJ,BRo/fi is the 
gyromagnetic ratio for electrons and is equal to 
2.21 X 10^ (rad.m).(A.s)“^, rb is the Bohr mag¬ 
neton, My = hqMsQ, and 1/At is proportional 
to the attempt frequency of the thermal field, At 
is the simulation time-step used, and the quantity 
G(o 7 )(t) is a Gaussian distribution with zero mean 

and unit variancd^SJ 

The thermal torque can be written aJ^ 


TTH(i) = My nm(t) X h(t) = Pgit) e^ — T’</,(t) eg 

( 10 ) 

where 


Pgit) = My[hxit) cosOit) cospit) 

+hyit) cos6it)sin(pit) — hzit) sinOit)], 

( 11 ) 

P(f,it) = My[hyit) cospit) — hxit) sinpit)]. (12) 

The magnetization dynamics under the action 
of the torques TE(t) and TTH(i) is described by the 
stochastic Landau-Lifshitz-Gilbert (LLG) equation 
as follows. 


dnm(t) ( dnm(t)\ 

—t — “ 

“-I/jITeW+TthWI- (13) 


After solving the above equation analytically, 
we get the following coupled equations of magneti¬ 
zation dynamics for 0 (t) and 


(1 -k a^) = ^[BQei<i>it))sineit) 

— 2aBi(pit),(jit))sin6it)cos9it) 

+ {aPgit) + P^m, (14) 


(l + «2) 


dpit) I 7 I 


[aBoeiPit)) 


dt My 
+ 2BiPit),ait))cos9it) 

- [sm6»(t)]“^ iPeit) - aP^it))] 

isin9 7 ^ 0). 


( 15 ) 
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We need to solve the above two coupled equations 
numerically to track the trajectory of magnetization 
over time, in the presence of thermal fluctuations. 

When sin0 = 0 (0 = 0° or 0 = 180°), i.e. when 
the magnetization direction is exactly along the easy 
axis, the torque on the magnetization vector given 
by Eq. ([6|) becomes zero. That is why only ther¬ 
mal fluctuations can budge the magnetization vec¬ 
tor from the easy axis. Consider the situation when 
9 = 180°. From Eqs. (flTl) and (fTKI) . we gell^^ 


4>{t) = tan 


/ ahy{t) -\- hx{t) \ 
\hy{t) - ahx{t)) ’ 


(16) 


de{t) 

dt 


+ hlit)) 


a/ {hy{t) - ah^{t)Y + {ahy{t) + K{t)Y 

(17) 

We can see from the Eq. (11711 clearly that thermal 
torque can deflect the magnetization from the easy 
axis since the time rate of change of 0[t) would 
be non-zero in the presence of the thermal fluctu¬ 
ations. Note that d9{t)/dt does not depend on the 
component of the random thermal field along the 
z-axis, i.e. hzit), which is a consequence of having 
2 -axis as the easy axis of the nanomagnet. However, 
once the magnetization direction is even slightly de¬ 
flected from the easy axis, all three components of 
the random thermal field along the x-, y-, and 2 - 
direction would come into play. 

When no stress is applied on the magnetostric- 
tive nanomagnet, magnetization would just fluctu¬ 
ate around an easy axis provided that the shape 
anisotropy energy barrier is enough high to pre¬ 
vent spontaneous reversal of magnetization from 
one state to another in a short period of time. We 
can solve the Eqs. (fT4]l and (fT5]l while setting Bstress 
= 0 to track the dynamics of magnetization due to 
thermal fluctuations. So this will yield the distri¬ 
bution of the magnetization vector’s initial orienta¬ 
tion when stress is turned on. Since the most prob¬ 
ably value of magnetization is along easy axis, the 
^-distribution is Boltzmann peaked at 0 = 0° or 
180°, while the (/>-distribution is Gaussian peaked 
at (j) = ±90° because these positions are minimum 

energy position^^fil 

Stress is ineffective when 6 is around 0° or 
180°, i.e. when magnetization is around the easy 
axis (mathematically, note that the expression in 
Eq. ([6]) is proportional to sin 9, which is zero when 
9 is equal to 0° or 180°). Hence, we would get a long 
tail in the switching delay distribution should mag¬ 
netization starts very near from easy axis. When we 


start out from 0 = 0°, 180°, we have to wait a while; 
thermal fluctuations may help here while getting 
started but random thermal kicks may also cause 
magnetization to traverse towards the opposite di¬ 
rection than the intended dirction of switching. 
Thus, switching trajectories initiating from near 
easy axis may be very slow. We do need to worry 
about the switching delay tail more than mean 
switching delay since the extent of tail will set the 
requirement of pulse width of stress for switching 
to take place with sufficiently high probability. 

In order to eliminate the long tail in the switch¬ 
ing delay distribution, we can apply a static bias 
field that will shift the peak of 9initiai distribution 
away from the easy axis, so that the most proba¬ 
ble starting orientation will no longer be the easy 
axi m. This field is applied along the out-of-plane 
hard axis (±a:-direction) and the potential energy 
due to the applied magnetic field can be expressed 
as 

Emag{9{t), cl){t)) = -MyH sin9{t) coscl){t), (18) 

where H is the magnitude of magnetic field. The 
torque generated due to this field is 

Tmit) = -nm(i) X VEmag{9{t),(t){t)). (19) 

We assume that a permanent magnetic sheet will 
be employed to produce the bias field and thus will 
not require any additional energy dissipation to be 
generated. The presence of this field will modify 
Eqs. ()14[) and ()15[) for magnetization dynamics to^^ 


d9{t) |7| 


dt My ' 

[BQe{(t>{l))sin9{t) — 2aB{(f){t), a{t))sin9{t)cos9{t) 
± aMvH cos9{t) cos4){t) — MyEd sin4’{t) 

+ + P4>{t))], ( 20 ) 


(1±«^) 


d(t){t) |7| 


dt My 
[aHoe(</>(t)) ± 2B{(j){t), a{t))cos9{t) 

— [sin9(t)]~^ [MyH cos9{t) cos<p{t) ± aMyH sin<p{t)) 
— [sin9{t)]~^ {PQ{t) — aP^{t))] {sin9^tS). (21) 


Note that the bias field makes the potential en¬ 
ergy profile of the nanomagnet asymmetric in (j)- 
space and the energy minimum gets shifted from 
= ±90° (the plane of the nanomagnet) to 


= cos 


-1 


H 


Ms{Nd_xX - Nd_yy)j 


( 22 ) 
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However, the potential profile will remain symmet¬ 
ric in 0-space, with 0 = 0° and 0 = 180° remaining 
as the minimum energy locations. A bias magnetic 
field of flux density 40 mT applied perpendicular to 
the plane of the magnet would make (j}min — ±87° 
deflecting the magnetization vector ~3° from the 
magnet’s plane. Application of the bias magnetic 
field will also affect the in-plane shape anisotropy 
energy barrier Ey as it gets reduced from 44 kT to 
36 kT at room temperature. 

We consider both the energy dissipated inter¬ 
nally in the nanomagnet due to Gilbert damping 
(termed as Ey) and the energy dissipated in the 
switching circuit while applying voltage across the 
multiferroic structure generating stress on the nano¬ 
magnet (termed as dissipation, where C and 

V denote the capacitance of the piezoelectric layer 
and the applied voltage, respectively). If the voltage 
is turned on or off abruptly then the energy dissi¬ 
pated during either turn on or turn off is (1/2)CI7^, 
however, if the ramp rate is finite, the energy dissi¬ 
pated can be significantly reducecP3. The internal 
energy dissipation Ey,, is given by the expression 
fo Pdit)dt, where r is the switching delay and Pd{t) 
is the power dissipated during switching given ai^^ 


^ (1+ l^)Mv ITEW+TMWy (23) 

We sum up the power Pd{t) dissipated during the 
entire switching period to get the corresponding en¬ 
ergy dissipation Ed and add that to the dis¬ 

sipation in the switching circuit to find the total 
dissipation Etotai- There is no net dissipation due 
to random thermal torque, however, it affects Ed 
since it raises the critical stress needed to switch 
with ~100% probability and it also affects the stress 
needed to switch with a given probability. 


2.2. Array of Multiferroic Devies 


Here we will use the same model as derived for 
a single multiferroic device and see how unidirec¬ 
tional flow of signal is possible in a horizontal chain 
of multiferroic devices^ 122] ^jg^ng dipole cou- 
nanomagnets and Bennett clocking 

The dipole coupling between two magnetic mo¬ 
ments Ml and M2 separated by a distance vector 
R can be expressed ajl 


pling between 
mechanisnJ2SI, 


P dipole — 


1 


47r/ro7?^ 


(M 1 .M 2 ) - ^(Mi.R)(M 2 .R) 

(24) 


1st 


2nd 


Nanomagnet Nanomagnet 


Magnetization 
'0 _ Direction 


0 


Magnet's plane: cf) = ± 90° 


Fig. 2. Dipole coupling between two magnetic moments. 
(Reprinted with permission from Ref. 21. Copyright 2012, 
Kuntal Roy.) 


In standard spherical coordinate system (see 
Fig. m, the expression of dipole coupling can be 
formulated as 

Edipole = [COS01COS02 

+ sin0isin62{cos4>icos4>2 — 2sin(^isin(f)2)] 

(25) 

where |Mi| = IM 2 I = roMsQ, Q is the volume of 
the nanomagnets. Mg is the saturation magnetiza¬ 
tion, and R = i?ey. 

Note that dipole coupling is bi-directional, i.e. 
Edipole — Edipole,! — Edipole, 2 ' Because of the dipole 
coupling between the magnetizations of the nano¬ 
magnets, the potential profiles of both the nano¬ 
magnets are tilted and the ground state of the mag¬ 
netizations are antiferromagnetically coupled as de¬ 
picted in the Fig. [2j 

If we somehow change the magnetization direc¬ 
tion of one nanomagnet, the magnetization of the 
other nanomagnet would not automatically change 
its direction to assume an antiferromagnetic order. 
It is because of the reason that there is a barrier 
separating two magnetization states. It’s true that 
antiferromagnetic order is the ground state, how¬ 
ever, during operation of devices, we must remove 
the barrier and then again restore it to make sure 
that antiferromagnetic order is maintained. Magne¬ 
tization may come to antiferromagnetic order after 
a very long time depending on the barrier height 
but the operation of devices cannot be dependent 
on that. 

In general, we need to propagate a logic bit uni- 
directionally along a chain of nanomagnets. It re¬ 
quires a clock signal to periodically reset the mag¬ 
netization direction of each nanomagnet. If a global 
magnetic field is utilized for such a purpose, it would 
not allow pipelining of data, and magnetization of 
every nanomagnet must be maintained along hard 
axis until a bit propagates. It needs an energy min¬ 
ima along hard axis, which can be introduced by bi- 
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4) = ± 90° 

2nd Nanomagnet: Energy Landscape 



Stress on 
Magnets 
2 and 3 



Remove 
Stress from 
Magnets 2 



Tie resolved 




Fig. 3. Imposing the unidirectionality in time to propagate a logic bit through a chain of nanomagnets. The 2nd and 3rd 
nanomagnets are stressed to align their magnetizations along the hard axis and then stress is removed/reversed on the 2nd 
nanomagnet to relax its magnetization along the desired state. (Reprinted with permission from Ref. 21. Copyright 2012, 
Kuntal Roy.) 


axial anisotropjSni^ thermal fluctuations would 
produce a large bit error probabilitjSB. Using a lo¬ 
cal magnetic field eliminates the problems of us¬ 
ing global magnetic field, bnt it is difficult to main¬ 
tain a magnetic field locally within a dimension of 
~100 nm. Furthermore, generating magnetic field 
is highly energy consnming. We can use electric- 
field operated (since electric-field can be maintained 
locally) multiferroic devices to propagate signals 
in a chain of nanomagnetJ^Q] using so-called Ben¬ 
nett clocking mechanism, termed in the name of 
BennetfP^. Ref. [20] performed the steady-state anal¬ 
ysis, while Refs. [2T1[22| solved the magnetization dy¬ 
namics using the same model as for a single multi¬ 
ferroic device to show that the switching may take 
place in snb-nanoseconcP^ ^33 [ 151 which is crucial 
for building nanomagnetic logio^^. 

Fig. [3] depicts the issue (and also solntion) be¬ 
hind Bennett clocking in a chain of nanomagnet. 
First of all, it needs to be emphasized that dipole 
conpling is bi-directional. So the 2nd nanomagnet 


experiences dipole coupling effect from both of its 
neighbors, i.e. 1st and 3rd nanomagnets. Note that 
we are considering only nearest neighbor interac¬ 
tion, since dipole coupling reduces drastically with 
distance [see Eq. ((Mj)]- Thus, if the 1st nanomag¬ 
net is switched, the 2nd nanomagnet finds itself in 
a locked condition as the 1st nanomagnet is telling 
it to go up, while the 3rd nanomagnet is telling it 
to go down. Therefore, it remains on its previons 
position and thus the change in information on the 
1st nanomagnet cannot be propagated throngh the 
chain of nanomagnets. 

To prevent this lockjam, we need to impose the 
unidirectionality in time as shown in the Fig. [3l 
Both the 2nd and 3rd nanomagnets are stressed to 
get them aligned to their hard axes (note the third 
row in Fig. [3]) and then stress is removed/reversed 
on the 2nd nanomagnet (note the fourth row in 
Fig. [3]) to relax its magnetization towards the de¬ 
sired state. In this way, snbseqnently applying stress 
on the nanomagnets and then releasing/reversing 
the stress, we can propagate a logic bit unidiree- 
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tionally along a chain of nanomagnets. The slight 
deflection in the magnetization of the 4th nanomag¬ 
net in the third row of Fig. [3] is due to dipole cou¬ 
pling, while in the fourth row, the magnetization of 
4th nanomagnet is aligned along its hard axis be¬ 
cause of applied stress on it. A 3-phase clock would 
be sufficient to propagate a signal along the chain of 
nanomagnet. Note that we are explaining the oper¬ 
ation with two-dimensional in-plane potential land¬ 
scapes of the nanomagnets (assuming azimuthal 
angle 4> = ±90°), but solution of the full three- 
dimensional dynamics is necessary since the out-of- 
plane excursion of magnetization has immense influ¬ 
ence in shaping the magnetization dynamicd^Sl and 
reducing the switching delay by a couple of orders 
in magnitude to sub-nanosecond^ 

The first row in the Fig. [3]shows a chain of four 
nanomagnets and we intend to switch the 2nd nano¬ 
magnet successfully in its desired direction once the 
1st nanomagnet is switched as depicted in the sec¬ 
ond row of Fig. [3l We will use subscripts 1-4 to 
denote the parameters and metrics for the corre¬ 
sponding nanomagnets. The dipole coupling energy 
on the 2nd nanomagnet due to 1st and 3rd nano¬ 
magnets can be written following the similar pre¬ 
scription given in the Eq. 


jm 


r^U 2 2 r 

Edipoie,2 = Mg Q [cos92Cos6i + cos 6 » 2 Cos 6»3 

± sin6isin92{cos(j)icos(j)2 — 2sin(/)isin(/>2) 
± sin9^sin92{cos(j)scos(j)2 — 2sin(pssin4>2)]- 

(26) 


and 


1 dEdi. 


' dipole ,02 


dipole,2 


sin92 d(j)2 
=_ 

\sin9i{cos4)isin4)2 ± 2sin(j)icos(j)2) 

± sin93{cos(j)3sin(f>2 ± 2sin(p3Cos(l)2)]- 

(29) 

The torque acting on the 2nd nanomagnet due to 
shape and stress anisot ropy can be derived similarly 
following the Eq. ([6]) aJ^D 

TE,2(i) = -2B2{(l)2{t))sin92{t)cos92{t) 

- Boe 2 {4>2 (t) ) sin92 {t) 60 , (30) 

where 


-\-Nd—yySiTl (^2(^) ^d—Zz\ 

+(3/2)A,a2^^, (31a) 

Boe2{4>2{t)) = Y^s^i^d-xx “ Nd_yy) sin{2<j)2{t)) ■ 

(31b) 

After solving the Landau-Lifshitz-Gilbert 
(LEG) equation considering the dipole coupling 
term in a very similar way as done for a single 
multiferroic device, we get the coupled dynamics 
between the polar angle 9 2 an d azimuthal angle (/>2 
for the 2nd nanomagnet as^ 


The torque ac ting on the 2nd nanomagnet due to 
dipole coupling2U 


Tdipole,2(^) — Ilm(^) ^ ^Edipole,2 


dE. 


dipole,2 


± 


1 dEdi 


dipole,2 




892 ^ sin92 d4>2 

Edipole,(j)2^(l> T Edipole,92^9 ■> (2^^) 


where 


^ 8Edipole,2 

dipole, 4>2 — o 7 i 
OU2 

2 2 r 

= -— MgQ. [—sin92Cos9i — sin92Cos93 
47rit° 

± sin9icos92{cos4>icos4>2 — 2sin4>isin(j)2) 
± sin93Cos92{cos(f>3Cos(f>2 — 2sin(p3sin(j)2)], 

( 28 ) 


(l ± = -^[Boe2iMi))sin92it) 

— 2 aB 2 {(l) 2 {i))sin 62 {t)cos 62 {t) 
T(iipole ,62 (32) 

± 2B2{4’2{t))sin92{t)cos92{t) 

± OlTdipole ,62 ± Edipole,(p2\ 

{sin9 / 0). (33) 

Note that in a very similar way the equations 
of dynamics for the other three nanomagnets can 
be derived. 

On energy dissipation, we have one more com¬ 
ponent contributing to the total energy apart from 
the shape anisotropy and stress anisotropy energy, 
which is the energy due to dipole coupling. While 
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calculating internal energy dissipation, the snm of 
the energy dissipations in all the fonr nanomagnets 
are considered but note that the dissipations in 1st 
and 4th nanomagnets are quite negligible since they 
don’t quite switch and dissipation in 2nd nanomag¬ 
net is around twice that of in 3rd nanomagnet since 
2nd nanomagnet switches a complete 180°, while 
the 3rd nanomagnet switches only about 90°. The 
instantaneous power dissipation for the 2nd nano¬ 
magnet can be determined aJ^ 


Pd^t) 


I 

Weak Bias 
Field 


Weak Bias 
Field 

I 


«l7l 

(1 + a‘^)Mv 


|TE,2(i) + Tdipole,2 

(34) 



(a) 


i t 


0 

1 

A B 

Out 

0 0 

1 

0 1 

1 

1 0 

1 

1 1 

0 



(b) 


I t 

0 1 


A B 

Out 

0 0 

1 

0 1 

0 

1 0 

0 

1 1 

0 


Fig. 4. Schematic of universal logic gates employing Mag¬ 
netic Quantum Cellular Automata (MQCA) based architec¬ 
ture: (a) NAND gate, and (b) NOR gate. Note that a weak 
bias field in the specified direction is required to break the tie 
when the input bits are different. The bias held must be weak 
enough so that it does not interfere in the operation when 
the input bits are Os for NAND gate and Is for NOR gate. 
(Reprinted with permission from Ref. 21. Copyright 2012, 
Kuntal Roy.) 


Note that we have not considered thermal flu- 
cutations and also have not applied any out-of-plane 
bias field so the term TM(t) term as in Eq. ([23|) is 
absent here. The power dissipations are integrated 
throughout the switching period to get the energy 
dissipation due to Gilbert damping. We have also 
considered energy dissipation, which can be 

significantly brought down by decreasing the stress 
since, stress is proportional to voltage applied, while 
sacrificing switching delay a bit. Since we have con¬ 
sidered instantaneous ramp and stress is reversed 
during ramp-down phase, the ‘CE^’ energy dissipa¬ 
tion is simply 3CV^ for the 2nd nanomagnet. 

We have considered Bennett clocking in an 
antiferromagnetically conpled horizontal wire for 
demonstration of magnetization dynamics in an ar¬ 
ray of multiferroic devices, however, a similar anal¬ 
ysis is possible in the context of a ferromagnetically 
coupled vertical wire. We have not incorporated 
ramp rate effect or thermal fluctuations, which one 
needs to consider and analyze further. Universal 
logic gates (e.g. NAND and NOR gates) can also 
be constructed and analyzed using the very same 
model that includes dipole coupling. Fig. |4] depicts 
such possibilities. In general, magnetizations of an 
array of nanomagnets can be manipulated to im¬ 
plement computing in MQCA (Magnetic Quantum 
Cellular Automata) based architecture^. Although 
such architecture has complexity of clocking each 
nanomagnet in the array, this is a regular structure 
and circuits based on such structure can be designed 
systematically. Anyway, unconventional design of 
logic gates and building blocks for large-scale cir¬ 
cuits using multiferroic composites can possibly be 
worked out too. Such designs may incur less com¬ 
plexity and possess better performance metrics than 
that of Bennett clocking mechanism. 


3. Simulation Results and 
Discussions 

In this Section, we review the simulation results for 
both single multiferroic deviceJ^ and an array of 
multiferroic deviceJ^H 122] xhe performance metrics 
switching delay and energy dissipation are deter¬ 
mined and trade-off between them is presented, i.e. 
if we want to make the switching faster, it would 
cost higher energy dissipation. Also, we determine 
the distributions of switching delay and energy dis¬ 
sipation, and number of successful switching events 
in the presence of room-temperature thermal fluc¬ 
tuations for a single multiferroic device. 
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3.1. Single Multiferroic Device 

We consider the magnetostrictive layer to be made 
of polycrystalline Terfenol-D, nickel, or coballP. 
Terfenol-D has 30 times higher magnetostriction co¬ 
efficient in magnitude and it has the following mate¬ 
rial properties - Young’s modulus (Y): 8x10^*^ Pa, 
magnetostrictive coefficient ((3/2)As): -|-90xl0“®, 
saturation magnetization (Mg): 8x10® A/m, and 
Gilbert’s damping constant (a): 0.1 (Refs. HU HU 
H5l I46p . For the piezoelectric layer, we have con¬ 
sidered lead-zirconate-titanate (PZT) having a di¬ 
electric constant of 1000. The maximum strain that 
can be generated in the PZT layer is 500 ppnPS SHI 
and it would require a voltage of 66.7 mV because 
(i3i=1.8xl0“^° m/V for PZT®. The PZT layer is 
assumed to be four times thicker than the mag¬ 
netostrictive layer so that any strain generated in 
it is transferred almost completely to the magne¬ 
tostrictive layeJ^QJ 0 gg corresponding stress 
in Terfenol-D is the product of the generated strain 
(500 X 10“®) and the Young’s modulus (8x10^® Pa). 
Hence, 40 MPa is the maximum stress that can 
be generated in the Terfenol-D nanomagnet. The 
strain-voltage relationship in PZT is actually su- 
perlinear since d^i increases with electric field®. 
Hence, the voltage needed to produce 500 ppm 
strain in the Terfenol-D layer will be less than 66.7 
mV and the energy dissipation would be a bit over¬ 
estimated too. 

We first review the results for low stress lev¬ 
els leading to slow switching speed (10-100 ns) and 
low energy dissipationP. Fig. [5] (taken from Ref. [9]) 
shows the energy dissipated in the switching cir¬ 
cuit (CV^) and the total energy dissipated (Etotai) 
as functions of delay for three different materials 
(Terfenol-D, nickel, and cobalt) used as the mag¬ 
netostrictive layer in the multiferroic nanomagnet. 
We solve magnetization dynamics to calculate the 
switching delay r and also energy dissipation (‘CV^’ 
dissipation and internal one Ey) for a given stress a, 
and then we plot the switching delays and energy 
dissipations for different stress values. Terfenol-D 
incurs much less energy dissipation than the other 
two materials because it has much higher magne¬ 
tostriction coefficient requiring a less stress level to 
generate a certain stress anisotropy. For Terfenol- 
D, the stress required to switch in 100 ns is 1.92 
MPa and that required to switch in 10 ns is 2.7 
MPa. Note that for a stress of 1.92 MPa, the stress 
anisotropy energy Bstress is 32.7 kT while for 2.7 
MPa, it is 46.2 kT. Since the stress level is assumed 


to be applied instantaneously, as expected, the en¬ 
ergy dissipation numbers are larger than the shape 
anisotropy barrier of ~32 kT. A larger excess en¬ 
ergy is needed to switch faster signifying the delay- 
energy trade-off. The energy dissipated and lost as 
heat in the switching circuit (CV^) is only 12 kT for 
a delay of 100 ns, while that is 23.7 kT for a delay 
of 10 ns. The total energy dissipated is 45 kT for 
switching delay of 100 ns and 70 kT for switching 
delay of 10 ns. 



Fig. 5. Energy dissipated in the switching circuit (CF^) and 
the total energy dissipated [Etotal) functions of delay for 
three different materials used as the magnetostrictive layer 
in the multiferroic nanomagnet. (Reprinted with permission 
from Ref. 9. Copyright 2011, AIP Publishing LLC.) 

With a nanomagnet density of 10^® cm“^ in a 
memory or logic chip, and if we consider 10% of the 
nanomagnets switch at any given time (10% activity 
level), the dissipated power density would have been 
only 2 mW/cm^ to switch in 100 ns and 30 mW /cw? 
to switch in 10 ns. Such extremely low power and 
yet high density magnetic logic and memory sys¬ 
tems can be powered by existing energy harvesting 
systemd® EB [52] [53| harvest energy from the 
environment without the need for an external bat¬ 
tery. These processors are deemed to be suitable for 
implantable medical devices, e.g. those implanted in 
a patient’s brain that monitor brain signals to warn 
of impending epileptic seizures. They can run on 
energy harvested from the patient’s body motion. 
For such applications, 10-100 ns switching delay is 
adequate. 

We now review multiferroic devices for higher 
stress levels and fast sub-nanosecond switching 
speed, which is particularly important for logic and 
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computing purposeJ^^. We performed simulations 
in the presence of room-temperature thermal fluc¬ 
tuations and we consider only Terfenol-D as magne- 
tostrictive material since it has much higher magne- 
tostrictive coefficient (than nickel and cobalt) thus 
being fruitful in resisting the adverse effects of ther¬ 
mal fluctuations. 

Fig. [6] plots the distributions of initial angles 
^initial and (pinitiai in the presence of thermal fluc¬ 
tuations and a bias magnetic field applied along the 
out-of-plane direction (-|-x-axis). The bias field has 
shifted the peak of BiniUai exactly from the easy axis 
{6 = 180°) as shown in Fig. [6Ka). The (piniUai dis¬ 
tribution (see Fig. EKb)) has two peaks and resides 
mostly within the interval [-90°,-|-90°] since the bias 
magnetic field is applied in the -l-x-direction. Be¬ 
cause the magnetization vector starts out from near 
the south pole {9 ~ 180°) when stress is turned on, 
the effective torque on the magnetization [ M X H, 
where M is the magnetization and H is the effec¬ 
tive field] due to the -l-x-directed magnetic field is 
such that the magnetization prefers the (^-quadrant 
(0°,90°) over the (^-quadrant (270°,360°), which is 
the reason for the asymmetry in the two distribu¬ 
tions of (pinitiai- Consequently, when the magneti¬ 
zation vector starts out from 9 ~ 180°, the initial 
azimuthal angle (pinitiai is more likely to be in the 
quadrant (0°,90°) than in the quadrant (270°,360°). 

We assume that when a compressive stress 
is applied to initiate switching (since Terfenol-D 
has positive magnetostrictive coefficient), the mag¬ 
netization vector starts out from near the south 
pole {9 ~ 180°) with a certain {9initial,(pinitiai) 
picked from the initial angle distributions at room- 
temperature. Stress is ramped up linearly and kept 
constant until the magnetization reaches x-y plane 
[9 = 90°). Then stress is ramped down at the same 
rate at which it was ramped up, and reversed in 
magnitude to aid switching. The magnetization dy¬ 
namics ensures that 9 continues to rotate towards 
0°. When 9 becomes < 5°, switching is deemed 
to have completed. A moderately large number 
(10,000) of simulations, with their corresponding 
{0initial,(pinitiai) picked from the initial angle distri¬ 
butions, are performed for each value of stress and 
ramp duration to generate the simulation results in 
the presence of thermal fluctuations. 

Fig. [7] shows the switching probability as a func¬ 
tion of stress levels (10-30 MPa) and as well as volt¬ 
age applied across the piezoelectric layer for differ¬ 
ent ramp durations (60 ps, 90 ps, 120 ps)l^ ^ at 


room temperature (300 K). The minimum stress 
needed to switch the magnetization with ~100% 
probability at 300 K is ~14 MPa for 60 ps ramp 
duration and ~17 MPa for 90 ps ramp duration. 
For higher ramp duration of 120 ps, the curve is 
non-monotonic, which we will discuss later. 


Distribution of e under thermai fluctuations 



(a) 

Distribution of <> under thermal fluctuations 



^ (degree) 
(b) 


Fig. 6. Distribution of polar angle Oinitial and azimuthal 
angle (pinitiai hue to thermal fluctuations at room tempera¬ 
ture (300 K) when a magnetic field of flux density 40 mT is 
applied along the out-of-plane hard axis (-|-a:-direction). (a) 
Distribution of polar angle Oinitial at room temperature (300 
K). The mean of the distribution is 173.7°, and the most 
likely value is 175.8°. (b) Distribution of the azimuthal an¬ 
gle (pinitiai due to thermal fluctuations at room temperature 
(300 K). There are two distributions with peaks centered at 
'^65° and ~295°. (Reprinted with permission from Ref. 19. 
Copyright 2012, AIP Publishing LLC.) 

At low stress levels (10-20 MPa), the switching 
probability increases with stress, regardless of the 
ramp rate. This happens because a higher stress can 
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more effectively counter the adverse effects of ther¬ 
mal fluctuations to facilitate switching, and hence 
increases the success rate of switching. This feature 
is independent of the ramp rate for lower stress lev¬ 
els. However, for higher stress levels accompanied 
by a higher ramp duration, the switching dynamics 
is complex, which we will explain onwards. 

Voltage (mV) 

16.7 25 33.3 41.7 50 



Fig. 7. Percentage of successful switching events among the 
simulated switching trajectories (or the switching probabil¬ 
ity) at room temperature in a Terfenol-D/PZT multiferroic 
nanomagnet versus (lower axis) stress (10-30 MPa) and (up¬ 
per axis) voltage applied across the piezoelectric layer, for 
different ramp durations (60 ps, 90 ps, 120 ps). The stress 
at which switching becomes ~100% successful increases with 
ramp duration. For large ramp duration (120 ps) or slow 
ramp rate, >^100% switching probability is unachievable. 
(Reprinted with permission from Ref. 19. Copyright 2012, 
AIP Publishing LLC.) 

We will now describe the significance of ramp- 
rate on success probability. When stress (compres¬ 
sive stress for Terfenol-D) is made active, mag¬ 
netization traverses from 6 ~ 180° towards x-y 
plane (9 = 90°). Once the magnetization vector 
crosses the x-y plane, the stress needs to be with¬ 
drawn as soon as possible. This is because the stress 
forces the energy minimum to remain at 9 = 90°, 
which will make the magnetization linger around 
9 = 90° instead of rotating towards the desired 
direction 0 ~ 0°. This is why stress must be re¬ 
moved or reversed immediately upon crossing the 
x-y plane so that the energy minimum moves to 
0 = 0°, 180°. Then the magnetization vector rotates 
towards 0 = 0° rather than 0 = 180° due to in¬ 
herent switching dynamicd^Sl. If the removal rate is 
fast, then the success probability is high since the 
harmful stress does not stay active long enough to 
cause significant backtracking of the magnetization 


vector towards 0 = 90°. However, if the ramp rate is 
too slow, then significant backtracking can possibly 
occur whereupon the magnetization vector may re¬ 
turn to the x-y plane. Then thermal torque decides 
the fate whether magnetization backtracks towards 
0 ~ 180°, causing switching failure or switches 
successfully towards 0 ~ 0°; in general there is 
50% switching probability then. Hence the switch¬ 
ing probability drops with decreasing ramp rate. 

Voltage (mV) 


16.7 25 33.3 41.7 50 



Fig. 8. The thermal mean of the switching delay (at 300 K) 
versus (lower axis) stress (10-30 MPa) and (upper axis) volt¬ 
age applied across the piezoelectric layer, for different ramp 
durations (60 ps, 90 ps, 120 ps). Switching may fail at low 
stress levels and also at high stress levels for long ramp du¬ 
rations. Failed attempts are excluded when computing the 
mean. (Reprinted with permission from Ref. 19. Copyright 
2012, AIP Publishing LLC.) 

A similar explanation is also applicable for the 
non-monotonic stress dependence of the switching 
probability when the ramp rate is slow (ramp du¬ 
ration of 120 ps). When 0 is in the quadrant [180°, 
90°], a higher stress is helpful since it provides a 
larger torque to move towards the x-y plane, but 
when 0 is in the quadrant [90°, 0°], a higher stress 
is harmful since it increases the chance of backtrack¬ 
ing, particularly when the ramp-down rate is slow. 
These two counteracting effects are the reason for 
the non-monotonic dependence of the success prob¬ 
ability on stress in the case of the slowest ramp 
rate. At higher stress levels accompanied by a slow 
ramp rate, it causes a significant amount of back¬ 
tracking causing the switching probability to drop 
fast. For slow ramp rate (ramp duration of 120 ps), 
we have not observed 100% switching probability at 
any stress for the 10,000 simulations performed. 

Fig. [8| plots the thermal mean (averaged from 
10,000 simulations) switching delay versus stress for 
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different ramp durations. Specifically, only success¬ 
ful switching events are considered here since we do 
not have a value of switching delay for an unsuc¬ 
cessful switching event. For ramp durations of 60 
ps and 90 ps, the switching delay decreases with 
increasing stress since the torque, which rotates 
the magnetization, increases when stress increases. 
However, for 120 ps ramp duration, the dependence 
is non-monotonic, due to same reason causing the 
non-monotonicity in Fig. [71 For a certain stress, 
decreasing the ramp duration (or increasing the 
ramp rate) decreases the switching delay because 
the stress reaches its maximum value quicker and 
hence switches the magnetization faster. 

Voltage (mV) 


16.7 25 33.3 41.7 50 



Fig. 9. The standard deviations in switching delay versus 
(lower axis) stress (10-30 MPa) and (upper axis) voltage ap¬ 
plied across the piezoelectric layer for 60 ps ramp duration at 
300 K. We consider only the successful switching events in de¬ 
termining the standard deviations. The standard deviations 
in switching delay for other ramp durations are of similar 
magnitudes and show similar trends. (Reprinted with per¬ 
mission from Ref. 19. Copyright 2012, AIP Publishing LLC.) 

Fig. [9] plots the standard deviation in switching 
delay versus stress for 60 ps ramp duration. The re¬ 
sults for other ramp durations are similar and hence 
are not shown for brevity. At higher values of stress, 
the torque due to stress is stronger and dominates 
over the random thermal torque, which causes the 
spread in the switching delay. This decreases the 
standard deviation in switching delay with increas¬ 
ing stress and thus the switching delay distribution 
gets more peaked as we increase the stress. 

Fig. Uni plots the thermal mean of the total 
energy dissipated to switch the magnetization as 
a function of stress and voltage across the piezo¬ 
electric layer for different ramp durations. To un¬ 
derstand these curves, we need to consider first 


the trend of average power dissipation {Etotai/f), 
which increases with stress for a given ramp dura¬ 
tion and decreases with increasing ramp duration 
for a given stress. More stress corresponds to more 
dissipation and also more internal power dis¬ 
sipation because it results in a higher torque. Slower 
switching decreases the power dissipation since it 
makes the switching more adiabatic. Note that the 
switching delay curves show the opposite trend (see 
Fig. [8]) than that of average power dissipation. At 
a lower ramp rate (higher ramp duration), the av¬ 
erage power dissipation Etotai/T is always smaller 
than that of a higher ramp rate, but the switching 
delay does not decrease as fast as with higher val¬ 
ues of stress (in fact switching delay may increase 
for higher ramp duration), which is why the energy 
dissipation curves in Fig. 1101 exhibit the cross-overs. 

Voltage (mV) 


16.7 25 33.3 41.7 50 



Fig. 10. Thermal mean of the total energy dissipation ver- 
sns (lower axis) stress (10-30 MPa) and (upper axis) voltage 
across the piezoelectric layer for different ramp durations (60 
ps, 90 ps, 120 ps). Once again, failed switching attempts are 
exclnded when computing the mean. (Reprinted with permis¬ 
sion from Ref. 19. Copyright 2012, AIP Publishing LLC.) 

Fig. dU plots the ‘CV^^ energy dissipation in 
the switching circuitry versus stress and the volt¬ 
age applied across the PZT layer. Increasing stress 
requires increasing the voltage V, which is why 
the energy dissipation increases rapidly with 

stress. This dissipation however is a small fraction 
of the total energy dissipation (< 15%) since a 
very small voltage is required to switch the magne¬ 
tization of a multiferroic nanomagnet with stress. 
The dissipation decreases when the ramp 

duration increases because then the switching be¬ 
comes more ‘adiabatic’ and hence less dissipative. 
This component of the energy dissipation would 
have been several orders of magnitude higher had 
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we switched the magnetization with an external 
magnetic fielcP13 or spin-transfer torque induced 
by spin-polarized curreniP. This makes the stress- 
induced switching of magnetization promising for 
pursuing our future nanoelectronics. 


Voltage (mV) 

16.7 25 33.3 41.7 50 



Fig. 11. The ‘CV^' energy dissipation in the external circuit 
as a function of (lower axis) stress and (upper axis) voltage 
applied across the PZT layer for different ramp durations. 
The dependence on voltage is not exactly quadratic since the 
voltage is not applied abruptly, but instead ramped up grad¬ 
ually and linearly in time. (Reprinted with permission from 
Ref. 19. Copyright 2012, AIP Publishing LLC.) 

Fig. [12] plots the switching delay and energy 
distributions in the presence of room-temperature 
thermal fluctuations for 15 MPa stress and 60 ps 
ramp duration. The high-delay tail in Fig. fT^ a) is 
in general associated with those switching trajecto¬ 
ries that start close to9 = 180°. In such trajectories, 
the starting torque is vanishingly small as explained 
earlier, which makes the switching sluggish at the 
beginning. During this time, switching also becomes 
susceptible to backtracking because of thermal fluc¬ 
tuations, which may also increase the delay further. 
However, it may well happen that random thermal 
torque quite facilitates a switching trajectory even 
it gets started very close to the easy axis making the 
net switching fast. Note that there was not a single 
event where the delay exceeded 1 ns out of 10,000 
simulations of switching trajectories, showing that 
the probability of that happening is less than 0.01% 
and probably far less than 0.01%. We have analyzed 
the cause of such high-delay tail and have been able 
to reduce it by application of an out-of-plane bias 
field. The magnitude of bias field can be calibrated 
to reduce the extent of such tail further. Since the 
energy dissipation is the product of the power dissi¬ 
pation and the switching delay, similar behavior is 


found in Fig. fT^b). 

Distribution of switching delay under thermal fluctuations 



Distribution of energy dissipation under thermal fluctuations 



Fig. 12. Delay and energy distributions for 15 MPa applied 
stress and 60 ps ramp duration at room temperature (300 
K). (a) Distribution of the switching delay. The mean and 
standard deviation of the distribution are 0.44 ns and 83 
ps, respectively, (b) Distribution of energy dissipation. The 
mean and standard deviation of the distribution are 184 kT 
and 15.5 kT at room temperature, respectively. (Reprinted 
with permission from Ref. 19. Copyright 2012, AIP Publish¬ 
ing LLC.) 


3.2. Array of Multiferroic Devices 

Here we review the simulation results of unidi¬ 
rectional information propagation through an ar¬ 
ray of multiferroic devices using Bennett clocking 
mechanisirPH fSH consider thermal fluc¬ 

tuations, nonetheless in all the simulations, the ini¬ 
tial orientation of the magnetization vector is as¬ 
sumed as 6 = 175° and cj) = 90°. Stress is applied 
instantaneously and we solve Eqs. (f32l) and (f33|) for 
the 2nd nanomagnet (and similar equations for the 
other three nanomagnets) at each time step. Once 6 
becomes 90°, stress is reversed instantaneously and 
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we follow the magnetization vector in time until 6 
becomes < 5°. At that point, switching is deemed 
to have completed. 

Fig. [13] shows the magnetization dynamics for 
all the nanomagnets when propagating a signal uni- 
directionally in a chain of four Terfenol-D/PZT 
multiferroic nanomagnets using Bennett clocking 
mechanism for the case of 5.2 MPa stress. When 
5.2 MPa compressive stress is applied on the 2nd 
and 3rd nanomagnets and after their magnetiza¬ 
tions come to their hard axes, stress is reversed 
on the 2nd nanomagnet to relax its magnetization 
towards its desired state. In the Fig. [Tdl bi. note 
that the applied stress has deflected the azimuthal 
angle of magnetization </> of the 2nd nanomagnet 
in the quadrant (90°, 180°) while that for the 3rd 
nanomagnet is deflected in the quadrant (0°, 90°). 
These out-of-plane excursions aid the switching to 
be very fast in sub-nanosecond; the same physics 
as for a single multiferroic devic J2IHHI applies here 
too. Had we not considered such out-of-plane ex¬ 
cursion, the switching delay would have been a 
couple of magnitudes larger, which signifies that 
steady-state analysis without considering out-of- 
plane motioEpSI is neither qualitatively nor quan¬ 
titatively accurat J2I] [22]_ 

The stress of 5.2 MPa has only shifted the mag¬ 
netizations out-of-plane by 5° during 90° switching, 
i.e. when the magnetizations come near to their 
hard axes. Upon reversing the stress on the 2nd 
nanomagnet, its magnetization rotates out-of-plane 
more, but this time at the very end, it has come 
back to its plane because of the dipole coupling with 
the 3rd nanomagnet, which tries to align the 2nd 
nanomagnet’s magnetization with its own magneti¬ 
zation. Finally, note that the magnetizations of 1st 
and 4th nanomagnets remain quite unchanged be¬ 
cause no stress is applied on these two nanomagnets, 
the slight changes in the directions of the magneti¬ 
zations occurred because of dipole coupling effect 
with the neighboring nanomagnets. 

Fig. 1141 shows magnetization dynamics for Ben¬ 
nett clocking with 30 MPa stress, rather than 5.2 
MPa stress as we have presented earlier. With high 
stress, magnetizations have deflected out-of-plane 
more (~10°) than that for the lower stress while 
reaching 9 = 90° and the magnetization of 2nd 
nanomagnet has executed a precessional motion be¬ 
fore completing switching. Apparently, almost half 
of the time taken during switching, magnetization 
experiences such unfruitful motion. To reduce such 


unfruitful motion, one can just use a lower stress 
level. If we apply 10 MPa stress, rather than the 
high stress 30 MPa, and magnetization does not 
execute any precessional motion and the switching 
delay increases just a bit from 0.32 ns to 0.34 ns, 
while energy dissipation decreases 10 times due to a 
lower stress level. Thus it is possible to engineer the 
stress levels to achieve good performance metrics 
for application purposes. 


Terfenol-D, 5.20 MPa Stress 



(a) 

Terfenol-D, 5.20 MPa Stress 



(b) 

Fig. 13. Magnetization dynamics for Bennett clocking in a 
chain of four Terfenol-D/PZT multiferroic nanomagnets with 
stress 5.2 MPa and assuming instantaneous ramp: (a) polar 
angle 9 versus time, and (b) azimuthal angle cj> over time while 
switching occurs, i.e. during the time 02 changes from 175° 
to 5°. (Reprinted with permission from Ref. 21. Copyright 
2012, Kuntal Roy.) 

Similarly, we can analyze and simulate logic 
gates (see Fig. m and bigger circuits using multifer¬ 
roic composites. Since the 3-phase clocking circuitry 
for Bennett clocking may be complex for application 
purposes, any unconventional design of logic gates 
and building blocks can be worked out and simu- 
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lated too, however, based on the same technique 
presented here. 



(a) 


Terfenol-D, 30.00 MPa Stress 



(b) 

Fig. 14. Magnetization dynamics for Bennett clocking in a 
chain of four Terfenol-D/PZT multiferroic nanomagnets with 
stress 30 MPa and assuming instantaneous ramp: (a) polar 
angle 0 versus time, and (b) azimuthal angle (f) over time while 
switching occurs, i.e. during the time 02 changes from 175° 
to 5°. (Reprinted with permission from Ref. 21. Copyright 
2012, Kuntal Roy.) 

4. Summary and Outlook 

We have theoretically investigated electric-field in¬ 
duced switching of magnetization in a magnetostric- 
tive nanomagnet strain-coupled with a piezoelectric 
layer in a multiferroic composite structure. The per¬ 
formance metrics of switching like sub-nanosecond 
switching delay and ~1 attojoule energy dissipa¬ 
tion are very promising for application purposes, 
which can extend the lifeline of conventional elec¬ 
tronics. Also we have performed simulations for an 
array of such multiferroic devices to demonstrate 
that unidirectional signal propagation is possible in 


sub-nanosecond switching delay, which is particu¬ 
larly important for general-purpose computing. In 
the same way, logic gates and building blocks for cir¬ 
cuits using multiferroic straintronic devices can be 
realized. Interestingly, the analysis of switching dy¬ 
namics shows both qualitatively and quantitatively 
that steady-state analysis is highly inaccurate par¬ 
ticularly in the context of switching process and de¬ 
termining switching delay. The out-of-plane excur¬ 
sion of magnetization highly facilitates the switch¬ 
ing process and it can reduce the switching delay by 
a couple of orders in magnitude to sub-nanosecond. 
So multiferroic straintronic devices are intriguing in 
respect to both applied physics and basic physics of 
binary switching. Utilizing multiferroic composites 
for the purposes of room-temperature computing 
can be so energy-efficient that it can be powered 
solely by energy harvested from the environment. 
Therefore, they are ideal for medically implanted 
devices which draw energy solely from the patient’s 
body movements, or even energy radiated by 3G 
networks and television stations. Moreover, the ba¬ 
sic building blocks are non-volatile permitting in¬ 
stant turn-on computer. This is an unprecedented 
opportunity in ultra-low-energy computing and can 
perpetuate Moore’s law to beyond the year 2020. 
Successful experimental implementation will pave 
the way for our future nanoelectronics. 
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